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Abstract 

Nonlinear time series analysis aims at understanding the dynamics of stochastic or chaotic processes. In recent years, 
quite a few methods have been proposed to transform a single time series to a complex network so that the dynamics 
of the process can be understood by investigating the topological properties of the network. We study the topological 
properties of horizontal visibility graphs constructed from fractional Brownian motions with different Hurst index 
H 6 (0, 1). Special attention has been paid to the impact of Hurst index on the topological properties. It is found that 
the clustering coefficient C decreases when H increases. We also found that the mean length L of the shortest paths 
increases exponentially with H for fixed length of the original time series. In addition, L increases linearly with 
respect to when H is close to 1 and in a logarithmic form when H is close to 0. Although the occurrence of different 
motifs changes with H, the motif rank pattern remains unchanged for different H. Adopting the node-covering box- 
counting method, the horizontal visibility graphs are found to be fractals and the fractal dimension dg decreases with 
H. Furthermore, the Pearson coefficients of the networks are positive and the degree-degree correlations increase with 
the degree, which indicate that the horizontal visibility graphs are assortative. With the increase of H, the Pearson 
coefficient decreases first and then increases, in which the turning point is around H - 0.6. The presence of both 
fractality and assortativity in the horizontal visibility graphs converted from fractional Brownian motions is different 
from many cases where fractal networks are usually disassortative. 
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1. Introduction 



Nonlinear time series analysis aims at understanding the dynamics of the underlying process. In recent years, 
numerous methods have been proposed to transform a single time series to a complex network |[lll2l - In this way, the 
dynamics of the process can be partially understood through the investigation of the network's topological properties. 
This direction is promising because there are a large number of tools developed in the field of network sciences 
i^mist]. Different transforming algorithms result in different networks, such as cycle networks based on the local 
extrema and their distances in the phase space fllziSli «-tuple networks lloi fioll. spa ce state networks based on 
conformational fluctuations | IT ] , visibility graphs based on the visibility of nodes lll2l llSIl , nearest neighbor networks 



in the phase space |l4 |15L 11611 . segment correlation networks lllTl Il8ll . temporal graphs 11191 120(1 . and recurrence 
networks 1211 122[ 12311 . In a recent work, a method is proposed to convert an ensemble of sequences of symbols into 
a weighted directed network whose nodes are motifs, while the directed links and their weights are defined from 
statistically significant co-occurrences of two motifs in the same sequence ll24ll . This method can also be applied to 
time series analysis. 
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Among the aforementioned methods, the visibility algorithm has attracted most applications from diverse fields, 
including stock market indices 1 25 , 2^ , human stride intervals [[ztIi , occurrence of hurricanes in the United States 1281. 
foreign exchange rates |29], energy dissipation rates in three-dimensional fully developed turbulence Isoll . human 
heartbeat dynamics llsil 32 1, electroencephalogram series llssll . binary sequences ll34ll . interevent time of book loans 
1 3511 . and daily streamflow series ll36ll . The visibility algorithm transforms a time series {x/j/^i ^ into a visibility 
graph G - {V, E), where V - {v,j,=i,...,iv is the set of vertex with the vertex v, corresponding to the data point x, and 
E = {eij]i,j=\,...,N is the adjacent matrix of the visibility graph whose element e,j = 1 if the following geometrical 
criterion is fulfilled: 



i — n 



I- J 



Vnj/ < n < j. 



(1) 



and etj - otherwise. The degree distribution of the visibility graph contains information of the original time 
series. For random series extracted from a uniform distribution in [0, 1], the degree distribution has an exponential tail 



P{k) ~ e ' " 111211 . In contrast, for fractional Brownian motions (FBMs), the degree distributions have power-law tails 
P{k) ~ k^", whose exponent a decreases linearly with the Hurst index H 11271 125 1 . 

Recently, a variant of the visibility algorithm has been proposed, which transforms time series into horizontal 
visibility graphs (HVGs) lll3ll . We note that the horizontal visibility algorithm is a variant of the visibility algorithm. 
A graph is an HVG if and only if it is outerplanar and has a Hamilton path fyf]. It has been proven that the degree 
distribution of any HVG mapped from random series without temporal correlations has an exponential form P{k) - 
(3/4)g-^(3/2)^ which allows us to distinguish chaotic series from independent and identically distributed (i.i.d.) time 



series 01311 . Furthermore, numerical simulations show that an HVG mapped from chaotic or correlated time series has 
exponential degree distribution 

P(k) ~ e-^'^, (2) 



in which a chaotic process has A < ln(3/2) and a correlated stochastic series has A > ln(3/2), separated by the i.i.d. 
case with A = ln(3/2) 

In this work, we investigate the topological properties of HVGs mapped from fractional Brownian motions with 
different Hurst indexes. Special attention is paid to the influence of the Hurst index of the fractional Brownian motion 
on the topological properties of the associated HVG. Specifically, the degree distribution, the clustering coefficient, 
the mean length of the shortest paths, the motif distribution, the fractal nature, and the mixing behavior of the HVGs 
are numerically studies. The most striking result is that the HVGs possesses both fractal and assortative features, 
which is different from the usual conclusion that fractal networks are disassortative 1 3^ 40ll . 

The paper is organized as follows. Section|2]studies the basic properties of the horizontal visibility graphs ma ppe d 
from FBMs. Section[3]investigates the fractal nature of the HVGs using the node-covering box-counting method 114 ill . 
Section|4]explores the mixing pattern of the HVGs. Section|5]summarizes our findings. 



2. Basic topological properties of the horizontal visibility graphs 

2.1. Illustrative examples of the horizontal visibility algorithm and HVGs 

The horizontal visibility algorithm transforms a time series {x,}i=i,...,A' into a horizontal visibility graph G - {V, E), 
where V - {v,} is a set of nodes corresponding to data points {x,j and E - {e,j) is the adjacent matrix of graph G, 
whose element eij - 1 if every points x„ between x, and xj fulfill the following condition lll3ll : 

Xi,Xj > x„, 'in\i <n < j. (3) 

Otherwise, e,j = 0, stating that point x, is not horizontally visible to point xj and the corresponding vertices v, and 
Vj are not connected. An illustrative example is shown in Fig.[T] FigurefTla) plots a simple time series containing 16 
data points by vertical thick bars, where the height of each bar is equal to the value of the corresponding data point. 
According to the horizontal visibility algorithm, we link every two bars when all the bars between them are lower 
than both of them and a horizontal link is drawn between these two bars that are visible to each other in the horizontal 
direction. The resultant horizontal visibility graph is shown in Fig.fllb). 

Figure|2]illustrates three typical HVGs mapped from FBMs with different Hurst indexes H - 0.9, 0.5 and 0.1 and 
the corresponding adjacent matrices. With the decrease of H, the HVG contains more loops and exhibits more fine 
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Figure 1: (Color online) The horizontal visibility algorithm, (a) A time series with horizontal visibility links, (b) The horizontal visibility graph 
mapped from the time series in (a). 

Structures and the diameter becomes shorter. We can conjecture from Fig. |2|a) that the mean length of the shortest 
paths is proportional to the size of the HVG for large H. By contrast, the HVGs for small H may exhibit small-world 
nature. In addition, the HVG mapped from an FBM with smaller Hurst index is more homogenous. These HVGs are 
reminiscent of fractal patterns such as viscous fingers, cracks, and diffusion-limited aggregates, implying that these 
graphs might have fractal nature. It is evident that the topological properties of the HVGs are strongly influenced by 
the Hurst indexes of the FBMs, which will be investigated in detail through numerical simulations. Consistently, the 
adjacent matrix varies in respect to the Hurst index. 
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Figure 2: Illustrative examples of HVGs. (a-c) HVGs mapped from FBM series of length N = 1000 with Hurst indexes H = 0.9, 0.5 and 0.1. (d-f) 
Adjacent matrices of the HVGs in (a-c). 



2.2. Degree distribution 

Degree distribution P{k) is one of the most important characteristic properties of complex networks. The degree 
distributions of the HVGs of correlated and uncorrelated stochastic and chaotic time series are exponential 1113 
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In addition, the visibility graphs of uncorrelated time series and fractional Brownian motions have exponential degree 
distributions and power-law tails, respectively Il2ll27ll25ll . A natural conjecture is that the degree distributions of the 
HVGs transformed from fractional Brownian motions might have power-law tails. We study the problem numerically. 
Figure [3] shows the degree distributions for the Hurst indexes H - 0.1, 0.3, 0.5, 0.7 and 0.9. It is found that the degree 
distributions have exponential tails as in Eq. (|2]i and the parameter A increases with H. 
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Figure 3: Empirical degree distributions of tlie hoiizontal visibility graphs converted from FBM series with different Hurst indexes H = OA, 0.3, 
0.5, 0.7, and 0.9. Each FBM series contains 10^ data points. The distributions have exponential tails and the parameter A increases with H. The 
solid line is the linear fit for H = OA. 



2.3. Clustering coefficient 

We now study the dependence of the clustering coefficient C of the HVG on the Hurst index of the FBM. The 
clustering coefficient C of a network defined as follows: 

C = (4) 

where 

where is degree of node v. The average of local clustering coefficient C states the probability of two neighboring 
nodes of any node are neighbors too. In other words, the clustering coefficient describes the occurrence frequency 
of triangles in the network (01 ■ Figure |4] shows that C decreases with H and the decrease speed also increases with 
H. This observation is evident according to Fig. [2] which shows that there are more triangles in the HVG when H 
decreases from 0.9 to 0.1. For small Hurst indexes, the clustering coefficients are large. 

2.4. Mean length of the shortest paths 

The mean length of the shortest paths of a network can be calculated as follows: 

' i=\ i=i+\ 

where d{i, j) is the length of the shortest path between node ; and j. For the HVG mapped from i.i.d. time series, the 
mean length of the shortest paths can be approximately derived: 

L(N) ^ 2ln N + 2{y- I) + 0(1 /N), (7) 
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Figure 4: Dependence of the clustering coefficient C of the HVGs against the Hurst index H of the original FBMs. Each FBM time series contains 
5 X lO' data points. For each H, the simulation is repeated for 100 times and the averages of C are shown. The inset shows the variance of C. 
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where y is the Euler-Mascheroni constant lll3[l . However, it is a hard task to obtain an analytic expression of L for 
correlated series 113 811 . Similarly, it is also difficult to derive the expression of L for the HVGs converted from FBMs. 
We thus study this problem by numerical simulations. 

According to Fig. [2] when H is close to 1, the HVG looks like a chain, which implies that the mean length L might 
be proportional to the FBM length A^. This is verified by numerical simulations. When H is close to 0, we find that L 
is proportional to In A^, indicating that these HVGs exhibit small-world feature. Therefore, we have 



L(A?) 



N, 



H 
H 



1, 
0. 



(8) 



As an example. Fig. |3a) illustrates the logarithmic dependence of the mean length of the shortest paths with respect 
to the HVG size for FBM series with H - OA. 




Figure 5 : (a) Logarithmic dependence of the mean length of the shortest paths (L) as a function of the HVG size (N) for FBM series with H = OA. 
For each N, 20 repeated simulations are performed and the average values of L are shown. The solid line is the linear fit of L versus In N to the 
data, (b) Exponential dependence of the mean length of the shortest paths (L) with respect to the Hurst index H for fixed HVG size. Three sizes 
N = 1000, 2000 and 3000 are illustrated. The solid lines are the linear fits of In L versus H to the three data sets. 
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As Fig. |2]suggests qualitatively, the mean length of the shortest paths increases with the Hurst index. We attempt to 
determine the quantitative relationship between L and H based on numerical simulations. The L values for different H 
are determined for fixed graph size A^. Figure|5lb) plots L with respect to H in semi-logarithmic scales for - 1000, 
2000 and 3000. It is observed that, for fixed A^, L increases exponentially in regard to H: 

L{H) ~ e)^", (9) 

where is a function of and independent of H. Figure|5jb) indicates that /3 increases with A^. 



2.5. Motif distribution 

Complex networks can be classified into different universality classes based on their macroscopic global properties 



1 4211 . Moreover, complex networks with the same global statistics may exhibit different local structure properties. At 



the microscopic level, the building blocks of complex networks are subgraphs or motifs and the network motif patterns 



in time series after mapping into nearest-neighbor networks, which is able to distinguish different types of dynamics 
in periodic, chaotic and noisy processes based on the occurrence frequency patterns of network motifs 11411 . For 



of occurrence can be used to define superfamilies of networks 11431 14411 . The superfamily phenomenon is also observed 

F dyn a 
Dl. 

nonstationary time series, the classification of superfamilies is strongly related to the DFA exponent of the time series 
ifl^ . It is thus interesting to investigate the pattern of motif distributions in the HVGs mapped from FBMs with 
different Hurst indexes. 

Since the HVGs constructed are undirected, we consider motifs of size 4. There are six different motifs as shown 
in Fig.|6|a). Our simulations are performed for different H values. For each H, we generate 20 fractional Brownian 
motions and convert them into 20 HVGs. For each HVG, the occurrence numbers of the six motifs are determined, 
denoted as n{M) where M - A, B, C, D, E and F. The occurrence frequencies of the six motifs are calculated as 
follows 

P(M) = n(M) / ^ n{M) . (10) 

' M=A,...,F 

For each H, the average values of P(M) are obtained. Figure |6jb) shows the motif distributions P(M) for different H. 
We find that the motif distribution patterns are same for all the H values: 

P(A) > PiB) > P{C) > P{D) > P{E) = P{F) = 0. 



For a given //, the motif occurrence frequencies for repeated simulations exhibit negligible fluctuations. Different 
from the nearest-neighbor networks lfl6ll . the Hurst index is not capable of distinguishing superfamilies. Plots (c-f) of 
Fig.|6]illustrate the nontrivial dependence of P{M) on H for M - A, B,C and D. 

The most striking feature of Fig.|6]is that n(E) — n{F) - 0. In other words, the two motifs E and F do not occur 
in all the simulated HVGs. Indeed, n{F) - holds for all HVGs and n{E) - almost holds for all HVGs. We provide 
a proof below. Consider four arbitrary data points {;<;,, , jc,, , x,, , x,,, } of a given time series, where i\ < (2 < i^ < u. For 
simplicity, we rewrite them as {x\,X2, X3, X4}. The four data points are mapped to four nodes of the resultant HVG. 

We first explain why motif F does not occur In motif F, each node is connected to the other three nodes. If there 
are any data points between x\ and X4, all of them must no larger than xi, X2, X3 and X4. Otherwise, it is impossible to 
transform {xi, X2, X3, X4} into a complete graph as motif F. Since there is a link between xi and X3, the X2 bar should 
be lower than xi and X3. We have X2 < X3. On the other hand, since X2 and X4 are visible, we have X3 < X2. These two 
constraints contradict, which means that the two hnks (xi, X3) and (x2, X4) cannot occur simultaneously. 

The explanation of «(£) = is slightly complicated. Since each node is connected to two nodes but unconnected 
to one node, there are three possible cases: (1) e,y = 1 except that e\2 - and 634 - 0, (2) e,j - 1 except that ei3 - 
and 624 = 0, and (3) e,j = 1 except that eu - and 623 — 0. For cases (1) and (3), en — 1 and 624 - 1- According to 
Eq. (O, we have X3 > X2 from ei3 = 1 and X2 > X3 from 624 - 1- Therefore, these two cases are forbidden. In case 
(2), we have ei4 = 1, ei3 - and 624 — 0. Since en - 1, all data points between xi and X2 are less than x\ and X2. 
Similarly, since 623 = 1, all points between X2 and X3 are less than X2 and X3. The fact that = results in x\ < X2 
or X3 < X2. In addition, we have X[ > X2 according to ei4 = 1. We thus have X3 < X2. Similarly, when we consider 
el4 = 1 and ^24 - 0, we have X2 < X3. At last, we obtain that X2 = X3. Combining all the information, motif E occurs 
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Figure 6: (a) AH six network motifs of size 4, which are labeled A, B, C, D, E and F, respectively, (b) Motif distributions for different FBMs with 
different Hurst indexes, (c) Dependence of occun'ence frequency of motif A as a function of H. (d) Same as (c) but for motif B. (e) Same as (c) but 
for motif C. (f) Same as (c) but for motif D. 

if and only if xi, X4 > X2 - xy, and ei2 - £23 - 634 = I. The absence of motif E in Fig.|6lb) simply means that the 
condition is not easy to be fulfilled for FBMs. 



3. Fractal scaling 

It is well known that many complex networks are small worlds ll45ll . Interestingly, small-world networks may 
also possess self-similar topological structure |41, 4^. The fractal dimension of the network can be obtained by the 



conventional box-counting method based on node covering |47, |48, 49, 50] or edge covering iBTIl . If a network is 
fractal, the minimal number Nb of boxes needed to cover the network scales with respect to the box size Ig as 141 1 



NbUb) ~ I 



-do 



(11) 



where cIb is the fractal dimension of the network characterizing its self-similar structure 115211 . On the other hand, 
scale-free networks may also exhibit self-similar connectivity behaviors after coarse graining processes, in which the 
degree distribution is scale invariant I53ll54ll55[l . Note that fractal scaling and self-similar connectivity behavior are 
two different properties ll56ll . In this section, we investigate the fractal nature of the HVGs based on the node-covering 
method. 



3.1. Box-counting method based on node covering using the simulated annealing algorithm 

According to the definition of fractal dimensions, we need to determine the minimal number of boxes of size / b to 



cover the nodes of a network. We use the simulated annealing algorithm 157115811 to determinate the minimum number 
Nb of boxes of size Ib needed to cover all the nodes of a network. The simulated annealing algorithm was used in 
the edge-covering method for the determination of the fractal dimension 15111 . The simulated annealing algorithm is 
implemented as follows. 

Starting from a given partition with C boxes of sizes no larger than Ig, we consider three types of operations to 
transform the partition into a new one: (1) One node is moved from one box with at least two nodes in it to another 
box if the diameters of both new boxes do not exceed /b; (2) One node is moved out of one box with at least two nodes 
to form a new box consisting of one node; and (3) Two boxes merge to form a new box, a move which is allowed if 
the diameter of the resulting box is no larger than Ib. 

At each temperature T, we perform k\ operations of the first type to exchange nodes, ko operations of the second 
type to form new boxes, and operations of the third type to merge boxes. An operation is accepted with probability 
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p, where 



P = 



1, Ca < Ch 

eXp[-iCa-Cb)/Tl Ca>Cb 



(12) 



where Cb and C„ are the numbers of boxes before and after an operation. It is obvious that an operation is more likely 
to be accepted at high temperature. After k\ + k2+ possible operations are performed, the system is cooled down 
to a lower temperature T' = cT, where c is a constant less than and close to 1 (typically c = 0.995 ~ 0.999). The 
typical starting temperature is around T = 0.6. If the number of boxes do not change after the ki + k2 + k^ operations 
at certain temperature, we argue that the algorithm reaches the optimal coverage and cease the annealing process. 

3.2. Fractal dimension ds 

We synthesize FBMs of size 3 x lO"' for different Hurst indexes. Each FBM is converted into a HVG. For different 
box sizes Is, the minimal numbers A^^ of boxes are determined using the node-covering method through simulated 
annealing as described in Section lTTI It is found for all HVGs investigated that Nb scales as a power law with respect 
to /fi, as expressed in Eq. (fTTT i. Figure |7j a) illustrates four typical examples of the nice power-law dependence with 
H - 0.1, 0.3, 0.5, 0.7 and 0.9. The linear regression of InA^^ against In Ib for each Hurst index gives an estimate of 
-dB- For each Hurst index, 10 FBMs are generated and 10 values of dB are determined. The dependence of the fractal 
dimension c/^ averaged over 10 realizations with respect to the Hurst index H is depicted in Fig.|7jb). We observe that 
the fractal dimension decreases with increasing H, which is consistent with the HVG structures shown in Fig.|2l 
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Figure 7: Fractal analysis of the HVGs. (a) Log-log plot of box number Ng against box size Ig for different FBM series with Hurst indexes 
H = 0.1, 0.3, 0.5, 0.7, and 0.9, respectively. The solid lines are best fits to the data for different Hurst indexes, (b) Linear relation between the 
fractal dimension dg and the Hurst index H when H > 0.3. Ten repeated calculations for each H are performed and the results are the averages. 
The dashed Une is f/g = 2 - // for comparison. 

From Fig.Qb), there is roughly a linear relation between ds and H when the Hurst index is larger than 0.3, and a 
linear regression gives that 

c/b = 1.96- 1.05//. (13) 

The error bar of the coefficient of the linear term is 0.03 and the corresponding /j-value is zero according to the t-test, 
showing that the coefficient is significantly different from zero and the linearity is confirmed. When H is very close 
to 1, we have dB ~ 0.91. It means that the HVG has a chain-like structure. When H is close to 0, we have dB ~ 1.96, 
which is however less than the extrapolated empirical value. At a first glance, it is surprising that the fractal dimension 
exceeds 2 while the HVGs are planar graphs. This reflects the fact that the planar HVGs for small //'s have to ruffle 
in the three- or higher-dimensional space so that the distance between two connected nodes is a finite constant while 
the distance between two unconnected nodes is infinite. It is thus possible to observe fractal dimensions larger than 2. 
For large Hurst indexes, it is interesting to compare the data with the following equation 

dB«2-H. (14) 
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It means that the fractal dimensions of these HVGs are close to, but less than, the fractal dimensions of the associated 
FBMs. This finding is consistent with the chain-like structures of the HVGs with large Hurst indexes. 



4. Assortative mixing pattern 

Mixing patterns are important statistical characteristics of complex networks ifsgilHoll . Assortative mixing refers 
to the phenomenon that nodes in a network are more likely connected to nodes that are alike them in some way. The 
mixing pattern can be measured using the Pearson coefficient 



(15) 



where M is the total number of edges, / = 1 , . . . , M, and y, and A:, denote the degrees of the two vertices at the ends of 
edge /. A network is assortative if r > and disassortative if r < 0. 
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Figure 8: Mixing pattern of HVGs. (a) Pearson coefficients r of HVGs mapped from different FBM series with different Hurst indexes H. The 
length of each FBM series is 10* and we simulated 30 times for each Hurst index H. Each data point represents the average value over the 30 
repeated simulations. The inset shows the standard deviation o"(r) of the Pearson coefficient r for different FBM seiies with different Hurst indexes 
H. (b) Dependence of the average degree {k„ (H)) of the neatest neighbors of nodes of the same degree k with respect to k for different H values. 

We synthesize FBMs of size - IC* with different Hurst indexes and transform them into HVGs. For each H, 30 
FBMs are simulated and 30 Pearson coefficients are calculated. The dependence of the average Pearson coefficient 
on the Hurst index is illustrated in Fig. [HI a). We find that r > for all H values, which is statistically significant in 
one standard deviation. In other words, we have r-cr,- > statistically. Therefore, almost all the HVGs mapped from 
FBMs are assortatively mixed with a few exceptions. 

An alternative measure of mixing patterns is the degree correlation For a given network, the average degree 
{kn{H)) of the nearest neighbors of nodes having degree k is determined. If {kn{H)) increases with k, high-degree 
nodes associate preferentially with other high-degree nodes and the network is assortative. Figure [8tb) shows the 
dependence of {k„{H)) with respect to k for different Hurst indexes. When H < 0.6, {kn{H)) increases linearly with k 
when A: > 3, which means that high-degree nodes are more likely to have high-degree neighbors. In addition, the slope 
of curve increases with decreasing H, which is consistent with the observation in Fig. [HIa) that r decreases with H. 
For the curves with H > 0.6, {k„{H)) increases when k < A and remains roughly constant when k > A, which means 
that low-degree nodes are more likely to have low-degree neighbors. In other words, high-degree nodes with A; > 4 are 
homogenously mixed, while low-degree nodes are assortatively mixed. This explains the large standard deviations of 
the Pearson coefficients for large H shown in the inset of Fig.[8|a). 
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For independent and identically distributed random time series, the average value of the degree k(x) of a node 
is a monotonically increasing function of the datum height x joll . Similar phenomenon is observed for fractional 
Brownian motions with different functional form, especially for H < 0.6. Therefore, high-degree nodes are more 
visible to each other, leading to the assortative mixing pattern in HVGs. 



5. Summary 

In recent years, ideas and tools in complex network sciences have been utilized to perform nonlinear time series 
analysis. A time series is transformed into networks based on various algorithms so that the dynamics of the process 
are mapped into the topological properties of the associated networks. A recently proposed algorithm based on 
horizontal visibility of data points is investigated in this work using fractional Brownian motions. 

We have studies the topological properties of the horizontal visibility graphs constructed from FBMs with different 
Hurst index H e (0, 1). It is found that the degree distribution, the clustering coefficient and the mean length L of the 
shortest paths all depend on the Hurst index. Specifically, with the increase of the Hurst index H, the absolute slope A 
of the exponential degree distribution increases, the clustering coefficient C decreases, and the mean length L of the 
shortest paths increases exponentially. In contrast, although the occurrence of different motifs changes with H, the 
motif rank pattern remains unchanged for different H. 

We also performed fractal analysis using the node-covering box-counting method based on the simulated annealing 
algorithm. The horizontal visibility graphs are found to be fractals and the fractal dimension decreases with 
H. Furthermore, the HVGs are found to exhibit assortative mixing patterns because the Pearson coefficients of the 
networks are positive and the degree-degree correlations increase with the degree. The assortativity of HVGs means 
that high-degree nodes are more likely connected with other high-degree nodes, es pec ially for H < 0.6. This feature 



cannot be explained by the hub-hub repulsion mechanism of fractality in networks 114711 . Our investigations show that 
hub-hub attraction can also produce fractality. It seems that the HVGs mapped from FBMs have different growth 
mechanisms compared with many real networks in which hub-hub attraction leads to non-fractality and assortativity 
and hub-hub repulsion leads to fractality and disassortativity 10]. 
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